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We report a high-resolution dielectric study on a pyrene-based discotic liquid crystal (DLC) in the 
bulk state and confined in parallel tubular nanopores of monolithic silica and alumina membranes. 

The positive dielectric anisotropy of the DLC molecule at low frequencies (in the quasi-static case) 
allows us to explore the thermotropic collective orientational order. A face-on arrangement of the 
molecular discs on the pore walls and a corresponding radial arrangement of the molecules is found. 

In contrast to the bulk, the isotropic-to-columnar transition of the confined DLC is continuous, 
shifts with decreasing pore diameter to lower temperatures and exhibits a pronounced hysteresis 
between cooling and heating. These findings corroborate conclusions from previous neutron and X- 
ray scattering experiments as well as optical birefringence measurements. Our study also indicates 
that the relative simple dielectric technique presented here is a quite efficient method in order to 
study the thermotropic orientational order of DLC based nanocomposites. 


I. INTRODUCTION 

The self-organization of disc-like molecules can result 
in liquid crystalline phases with remarkable structural 
and dynamical properties (THThI . In particular, the elec¬ 
trical conductivity and the high photoconductivity typ¬ 
ical of columnar phases make discotic liquid crystals 
(DLCs) promising candidates for numerous technologi¬ 
cal applications, for example in electroluminescent and 
photovoltaic devices fl9l-[23j. 

The molecular building blocks of many DLCs are pol¬ 
yaromatic cores surrounded by aliphatic chains. The 
stacking of the polyaromatic cores results in an hexagonal 
columnar structure. The one-dimensional conductivity 
originates in charge transfer along the molecular columns 
of the 7T-7T systems, whereas the alkyl side-chains insu¬ 
late each column from its neighbors [24f | . The columnar 
structure, therefore, represents an array of isolated con¬ 
ducting molecular wires. In comparison to conventional 
crystalline or amorphous conductors, liquid crystals have 
a number of advantages, such as simple processability, 
self-healing and spontaneous alignment. 

DLCs can also easily be embedded into nanoporous 
media by melt infiltration (25j - l3l| . One of the promising 
ideas in this respect is the synthesis of nanowires by this 
technique as pioneered by Steinhart et al [32|. However, 
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its practical implementation reveals difficulties to provide 
a long-range molecular order, which is easily maintained 
in the bulk state [13]. Possible strategies to overcome 
this problem are modifications of host-guest interactions, 
i.e. chemical modification of the pore wall along with a 
proper choice of the guest DLC material and the pore 
diameter 30], or controlling the packing structure by ap¬ 
plying external fields, e.g. by magnetic fields[33]. 

However, similarly as the dynamical and structural be¬ 
havior of rod-like liquid crystalline systems (34l - [47l ]. the 
properties of DLCs has turned out to be particularly 
susceptible to nano confinement |27H32l I48I - |50| . Thus, 
new, simple experimental methods that are able to mon¬ 
itor the type of molecular ordering inside nanoporous 
media are strongly demanded. Here, we demonstrate 
that a high-resolution dielectric study can yield impor¬ 
tant insights into the thermotropic, collective orienta¬ 
tional order of a DLC and thus on confinement effects 
on its phase behavior. We compare those findings with 
results from experimentally more challenging and elab¬ 
orated techniques, i.e. temperature-dependent neutron 
and X-ray diffraction [|51] as well as optical birefringence 
experiments l~2?s. l~50]. 

The DLC used in our studies was the pyrene-1,3,6,8- 
tetracarboxylic rac-2-ethylhexyl ester (Py4CEH). It con¬ 
sists of a polyaromatic core surrounded by flexible 
aliphatic chains (52}. In the bulk state Py4CEH forms 
a hexagonal columnar discotic (CD) phase (see Fig. la) 
between the isotropic (I) phase, which exists above 365 K, 
h 'ttaH^Wd 1 flfhstic crystalline phase observed below 246 K 
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FIG. 1. (color online) Illustrations of the capacitor geometry used for the dielectric measurements in the bulk state (a) and 
for the nano-confined state (b). Reduced dielectric permittivity of Py4CEH vs. temperature during heating (circles) and 
subsequent cooling (squares) measured in the bulk state (c) and confined to mesoporous matrices of different pore radii of 
PAI 2 O 3 (d,e,f) and pSi 02 (g,h,i). The dashed lines are e,. (^-dependences extrapolated from the isotropic phase. Ae r indicate 
the excess permittivity. 



II. 


EXPERIMENTAL 


We employed monolithic membranes of alumina 
(PAI 2 O 3 ) and silica (pSi 02 ), which are permeated by ar¬ 
rays of parallel, tubular channels. The PAI 2 O 3 mem¬ 
branes (thickness h =100 /rm, Smart Membranes GmbH 
(Germany)) have been fabricated by electrochemical 
etching of pure aluminum in a sulphuric acid electrolyte. 
Their pore size has been determined by volumetric N 2 - 


sorption isotherms recorded at T = 77 K which gives the 
following average pore radii, R: 21.0±2.0 nrn (porosity 
P=24%), 15.5T1.5 nm (P =17%) and 9.3T0.7 nrn (P 
=16%) - see the top-view electron microgrograph of the 
membrane with 19 nm mean pore diameter in Fig. 2d. 
The pSi 02 membranes have been fabricated by oxida¬ 
tion (12 h at T=800 °C under standard atmosphere) 
of porous silicon obtained by electrochemical etching of 
boron doped silicon wafers (resistivity 0 . 01 - 0.02 ohm-cm) 
using an electrolyte mixture of HFiCtHrOH (2:3) and a 
DC current density of 11-13 mA/cm 2 |54h 56| . By varying 
the etching time and etching current, we have fabricated 
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pSi0 2 membranes with pore radii 7.0±0.5 nm (P = 52%, 
h =340 /im ), 5.6±0.4 nm (P = 35%, h =290 /*m) and 
4.6±0.3 nm (P = 30%, h =180 /im). 

Dielectric measurements have been performed with the 
Frequency Response Analyzer Solatron-1260A. For the 
bulk measurements a sample cell made of two parallel 
ITO-glass electrodes separated by a gap, h = 10 /im, 
was used. The DLC has been aligned homeotropically, 
as sketched in Fig. la. 

A schematics of the capacitor design used in the di¬ 
electric measurements of the nanocomposites is shown 
in Fig. lb. Gold electrodes have been evaporated onto 
the membranes. Before filling of the samples they were 
annealed for 15 minutes at ~ 400 K in order to remove 
adsorbed pore water. Then the mesoporous membranes 
were completely filled by spontaneous imbibition (infil¬ 
tration) in the isotropic phase at 370 K (57j |. 

Note that in our capacitor geometry the channels are 
all arranged parallel to the applied electric field. The 
sample forms a simple parallel electrical circuit. Accord¬ 
ingly, the effective complex permittivity of the composite 
e* is given by the permittivities of the host matrix and 
the guest DLC e* d weighted with the their volume frac¬ 
tions: e* = (1 — P)e^ n -\-Pe* d . The temperature-dependent 
dielectric measurements have been performed with cool¬ 
ing and heating rates of 0.07 K/min. 

The temperature changes of the host pAl 2 C >3 or pSi0 2 
matrices are weak and smooth in the entire tempera¬ 
ture range with negligible losses in the frequency window 
10 2 < / < 10 6 Hz. Thus, the molecular ordering inside 
the pores caused by the isotropic-to-columnar hexago¬ 
nal phase transition, which particularly influences the e* d 
permittivity, can be easily detected by analyzing the tem¬ 
perature dependence of the permittivity of the nanocom¬ 
posite, e*(r)[58j. 


III. RESULTS AND DISCUSSION 

For the Py4CEH based nanocomposites the dielectric 
losses at frequencies below 1 kHz are due to conductivity 
related processes. These processes are not in the focus 
of this study. At higher frequencies, particularly close to 
the upper limit of the frequency window of the frequency 
response analyzer (/ < 10 6 Hz) the dielectric losses of 
confined DLC Py4CEH are comparable to the experi¬ 
mental resolution used here [59} . Because of the Kramers- 
Kronig relationship between dielectric loss and real part 
of the permittivity, the real part, e'(/), can be considered 
frequency-independent. This implies that the dipolar re¬ 
laxation corresponds to the so-called quasi-static limit, 
r -C (27r/) -1 , where r is the corresponding relaxation 
time. This is in agreement with recent dielectric stud¬ 
ies on bulk Py4CEH [53, ,59| , where the Arrhenius plot 
(see Fig. 4 of Ref.[59j) at the temperature region rele¬ 
vant here (from 310-370 K) give relaxation frequencies, 
f r = in the range between 30 and 500 MHz. 

Our measurements have been performed at a constant 


frequency of 200 kHz. This choice has the advantages 
that it results in a quasi-static character of the dielectric 
measurements (/ -C f r ) as well as an optimized accuracy 
and resolution, which depending on the measured capac¬ 
itance changes from ~ 10 -4 to 10~ 3 pF (i.e. 0.003 to 
0.01% of the absolute values). As we will document be¬ 
low, this high resolution appears particularly crucial for 
our study, given the relatively weak changes of the dielec¬ 
tric permittivity associated with orientational ordering of 
the discotic molecules. 

The absolute magnitude of the permittivity of the nano 
composite scales with the porosity and is, itself, not im¬ 
portant for our analysis. Therefore, it is convenient to 
characterize the dielectric temperature behavior in the 
region of the phase transition by the reduced dielectric 
permittivities, e r (T) = e'(T)/e'(T dCD ) (bulk state) or 
e r (T) = e' (T)/ e'(Tj?) (confinement), expressed simply by 
the ratio of the measured capacitances C'(T)/C(T dCD ) 
or C'(T)/C(Tfc), respectively. Here, T dCD and rep¬ 
resent characteristic temperature points of the isotropic- 
to-columnar transition (discussed below) during heating 
in the bulk state and in confinement, respectively. 

In Fig. 1 we present the temperature dependences of 
the reduced dielectric permittivity, e r (T), measured dur¬ 
ing heating and cooling runs for discotic Py4CEH in the 
bulk state and embedded into mesoporous matrices of 
different types and channel sizes. 

Compared to the liquid crystalline phase, the bulk per¬ 
mittivity, e r , in the isotropic phase has only a weak tem¬ 
perature dependence. The isotropic-to-columnar tran¬ 
sition is accompanied by a sudden rise of the reduced 
permittivity with a characteristic kink in the e r (In¬ 
dependence at Tj CD , see Fig. lc. Subsequent cooling 
reveals a small temperature hysteresis (~ 0.7 K) of the 
reverse transformation. In the bulk state the homeotrop¬ 
ically aligned columnar phase exhibits a positive ex¬ 
cess permittivity, Ae r . This may appear surprising at 
first glance, if one considers that at optical frequen¬ 
cies Py4CEH is characterized by a negative dielectric 
anisotropy (£|| i00 < £_l,oo) which leads to a negative bire¬ 
fringence in the columnar discotic phase [281 l50j |. It can, 
however, be rationalized as follows: The dielectric prop¬ 
erties at optical frequencies are mainly determined by 
the electronic polarizability. Here the contribution of the 
polyaromatic cores is dominant, which is characterized 
by a large polarizability anisotropy. At low frequencies, 
the dipolar polarizability and the orientational polariza¬ 
tion of molecular dipoles provide additional contributions 
to the dielectric permittivity. Accordingly, in the quasi 
static limit the dielectric permittivities are determined by 
£||,st = £||,oo + Ae|| and £j_ iSt = £_l i00 + Ae ± , where Ae|| 
and Ae_l are the dipolar relaxation strengths. A positive 
dielectric anisotropy of Py4CEH DLC means that an in¬ 
equality Ae || — Ae_l > £j_,oo ~ £iioo holds. Interestingly, 
the polyaromatic core does not contribute to the molec¬ 
ular dipole moment. Thus, the dipolar polarizability is 
attributable to the polar ester groups and the aliphatic 
chains forming the intercolumnar network in the hexago- 
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FIG. 2. (color online) Molecular ordering types of DLCs in 
cylindrical nanochannels. (a) Parallel axial configuration, (b) 
Radial, face-on configuration with isotropic core, (c) Radial 
multidomain face-on configuration, (d) Top-view electronmi- 
crograph of an alumina membrane with approx. 19 nm mean 
pore diameter. The scale bar indicates a width of 100 nm. 


nal liquid crystalline structure. Their flexibility and abil¬ 
ity to order as a function of temperature explain, presum¬ 
ably, why the quasi-static dielectric permittivity does not 
saturate at lower temperatures. This is in contrast to the 
saturation of the optical birefringence [H, [50] and thus 
to the saturation of the collective orientational ordering 
of the polyaromatic cores. 

The temperature behavior of the reduced permittivity, 
e r (T), for Py4CEH embedded into nanopores substan¬ 
tially differs from the bulk behavior as one can deduce 
from Figs.ld-li. For all nanocomposites studied, the ex¬ 
cess permittivity, defined as the difference between the 
measured value, e r (T) and the one obtained by an ex¬ 
trapolation of the e r .(T)-dependence from the isotropic 
phase (see marked by Ae r in Fig. Id for explanation), 
is found to be negative. This means that the parallel 
axial columnar order (see sketch in Fig. 2a), which is 


particularly desirable for many nanotechnological appli¬ 
cations, is not realized. The negative Ae r is consistent 
with a molecular ordering for which discotic molecules 
are oriented preferably parallel to the long pore axis. The 
physico-chemical reason for the realization of such an or¬ 
dering is the face-on type anchoring enforced by alumina 
and silica surfaces. It can lead to several different config¬ 
urations, the most probable are the radial configuration 
with isotropic core (Fig. 2b) and the radial multidomain 
configuration (Fig. 2c) [28l. l5ll [60]. 

This observation is consistent with our X-ray and neu¬ 
tron diffraction study on Py4CEH in alumina and silica 
channels j5jJ. It indicated that the molecular ordering 
is either of radial columnar or of radial nematic type 
for small channel diameter, see the X-ray and neutron 
diffraction patterns in Fig. 3 and Fig. 5, respectively, and 
the corresponding discussion in Ref. [5l[ ■ 




R ' 1 [nm' 1 ] 


FIG. 3. (color online) (a): The characteristic temperatures 
Tft and Tjr determined from the temperature dependences of 
the dielectric permittivity (see Fig. 1) measured during heat¬ 
ing (open circles) and cooling (solid circles), respectively. The 
dasli-dotted line is a guide for the eye. (b): Tj? vs. the in¬ 
verse channel radius R as determined by dielectric measure¬ 
ments (open circles) and by optical birefringence data (open 
squares), respectively, in comparison with a 1/R-fit (dash-dot 
line) and with a function T]f = T^ CD — a\R~ 1 — 02 -fiX 2 and the 
fit-parameters, Tf CD = 365.2 K, ai=52.1 K-nm and 02=241.0 
K-nm 2 (solid line), see Ref.[2^] and discussion in the text. 

As is obvious from Fig. 1 the geometrical constraints 
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render the transition in the nanochannels gradual and 
hysteretic with characteristic lower and upper closure 
points (kinks) in the dielectric data at Tj: (cooling) and 
Tfc (heating). Just below these characteristic temper¬ 
atures the excess permittivity, A e r (T), exhibits strong 
variations with a saturation at lower temperatures. The 
larger the channel diameter, the steeper are the changes 
of A e r {T). The characteristic temperatures, T£ and Tfc, 
decrease with decreasing channel radius. 

Note that, as discussed extensively in Refs. [H, H3] 
and Tj} do not represent phase transition tempera¬ 
tures in a classical meaning. The radial arrangement 
illustrated in Fig. 2b results in curved molecular layers, 
corresponding splay deformations and excess free ener¬ 
gies upon collective molecular orientation. These excess 
energies are not homogeneous in pore space. Because 
of the increase in the curvature of the layers with de¬ 
creasing distance r from the channel centre, they decrease 
with 1/r. Hence, the transition temperature towards the 
orientationally ordered state depends on the radial posi¬ 
tion of the layer. It decreases with decreasing r. During 
cooling the transition towards the radially ordered state 
starts at the periphery, near the channel walls, where 
the excess free energy is the lowest. Therefore, the in¬ 
terfacial interaction and the channel roughness play an 
important role for the onset of the phase transition. Ra¬ 
dial homeotropic and axial planar ordering may simul¬ 
taneously occur, enforced by specifics of the interfacial 
and/or local geometrical conditions. Accordingly, first a 
(metastable) coexistence of both states appear and only 
at lower temperature the radial molecular ordering be¬ 
comes dominating. With these assumptions we could ex¬ 
plain in our previous optical polarimetry study that the 
kink position upon cooling varied from sample to sample 
and was sensitive to the surface coating 28]. 

By contrast, the kink position during heating exhib¬ 
ited a more systematic evolution as a function of chan¬ 
nel diameter and was independent of the chemistry of 
the channel surface. Presumably, there is always an en¬ 
tirely disordered or highly defect decorated core - see 
Fig. 2b. After all, the splay-deformation related excess 
energy diverges towards the pore center. This disordered 
core acts as a nucleation for the high-temperature phase. 
The phase front starts there and gradually propagates 
according to the r-dependent transition temperature to 
the periphery, while ’’melting” subsequent radially or¬ 
dered molecular layers. The kink position at heating 
corresponds therefore to the situation at which the phase 
front reaches the channel wall. It is defined mainly by a 
geometrical factor, i.e. the channel radius, and because 
of the omnipresent nucleation mechanism in the channel 
center it is a reliable measure to characterize the phase 
transition. For channels of radius R with ideally smooth 
walls its temperature can be derived from a Landau-de 
Gennes free energy functional which includes the pecu¬ 
liar influence of radial-dependent excess energy outlined 
, i.e. Tk = Tf rn — a\R~ x — a, 2 R~ 2 ■ Here 


above 

Tiq D = 365.2 K is the bulk transition temperature, the 


1 /R term is analogous to the Gibbs-Thomson prediction 
for a phase transition of first order in cylindrical confine¬ 
ment [6ll - [63i] and the 1 /R 2 scaling results from the splay 
deformation contribution. 

In excellent agreement with these considerations also 
the T-shifts of the upper closer points of the hysteresis 
in the dielectric data, Tf} are dominated by two types of 
scalings, see Fig. 3: An 1/R-scaling at large channel di¬ 
ameters, where the splay deformation energies are negli¬ 
gible, and an additional l/i? 2 -scaling for smaller channel 
diameters, which is characteristic of the excess energies 
related to splay distortions in the collective arrangement 
of the discotic molecules [Hj]. Moreover, the prefactors of 
the scaling laws deduced from the Independence of the 
closure point shifts agree with the ones derived in the 
birefringence study, i.e. a± = 52.1 K-nm and 02 = 241.0 
K-nm 2 . 


IV. SUMMARY 

In summary, we have reported a high-resolution di¬ 
electric study on the columnar DLC Py4CEH confined 
in tubular nanopores. In contrast to the optical bire¬ 
fringence, the dielectric permittivity of Py4CEH in the 
quasi-static limit is characterized by a positive dielec¬ 
tric anisotropy in the columnar discotic phase. This 
determines the dielectric behavior in the region of the 
isotropic-to-colunmar discotic transition both in the bulk 
and in the confined states. Whereas the homeotropic 
alignment in the bulk state is characterized by a posi¬ 
tive excess permittivity, molecular ordering of Py4CEH 
in the confined cylindrical geometry is accompanied by 
a negative value. It corresponds to ordered structures, 
where discotic molecules orient preferably parallel to the 
long channel axis (in a so-called face-on configuration). 
The R —scaling of characteristic phase transition temper¬ 
ature points, deduced from the dielectric data, are in very 
good agreement with recent neutron scattering and opti¬ 
cal birefringence studies on the same discotic mesogen. 

Molecular dynamics simulations on rod-like systems 
have successfully contributed to a more detailed under¬ 
standing of the isotropic-to-nematic transition in confine¬ 
ment [3a. Ell, ffijtiiifii] . In fact, a recent molecular dynam¬ 
ics study documented radial movements of phase fronts 
analogous to the ones inferred here, however for a rod-like 
system in nanochannels [tI|. Molecular dynamics simu¬ 
lations on confined and bulk DLC [60, |67| indicate, in 
agreement with the experiments presented here, that the 
phase behavior of the confined system strongly differs 
from its bulk counterpart: the bulk isotropic-to-columnar 
transition is replaced by a continuous ordering from a 
paranematic to a columnar phase. Moreover, a reorga¬ 
nization of the intercolumnar order was observed which 
reflects the competing effects of pore surface interaction 
and genuine hexagonal packing of the columns. 

This competition between surface anchoring and ten¬ 
dency to form rigid, hexagonal columnar phases has 
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turned out to be decisive for the formation of either ax¬ 
ial or radial or circular concentric arran gem ents of the 
columns in cylindrical nano confinement {2714321 148l - f50i | . 
For edge-on anchoring Zhang et. al. [29|, Ho[ found a tran¬ 
sition from circular concentric towards an axial colum¬ 
nar arrangement upon increase in columnar rigidity of 
a triphenylene based DLC. This result corroborates the 
observations of Duran et. aZ.[27[. They found for a 
DLC with a comparably large core (hexabenzocoronene), 
and thus presumably large column rigidity, confined in 
tubular alumina nanopores an axial arrangement of the 
columns in the channel center, despite a face-on anchor¬ 
ing at the pore walls. Moreover, below a certain pore di¬ 
ameter the transition towards a low-T crystalline phase, 
typical of the bulk of this system, is completely sup¬ 
pressed. 

Diffraction experiments and Molecular Dynamics sim¬ 
ulations on partial fillings [45], [68] may allow for com¬ 
plementary structural information on the propagation of 


the phase boundaries in channel space. This will shed 
light on the remarkably different pathways found here 
for freezing and melting of discotics. 
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